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The PHENIX experiment presents results from the RHIC 2005 run with polarized proton collisions 
at y/s = 200 GeV, for inclusive n° production at mid-rapidity. Unpolarized cross section results are 
given for transverse momenta px = 0.5 to 20 GeV/c, extending the range of published data to both 
lower and higher pr- The cross section is described well for pr < 1 GeV/c by an exponential in 
Pt, and, for pr > 2 GeV/c, by perturbative QCD. Double helicity asymmetries All are presented 
based on a factor of five improvement in uncertainties as compared to previously published results, 
due to both an improved beam polarization of 50%, and to higher integrated luminosity. These 
measurements are sensitive to the gluon polarization in the proton, and exclude maximal values for 
the gluon polarization. 

PACS numbers: 13.85.Ni,13.88.+e,21.10.Hw,25.40.Ep 



A principal goal of the spin program at the Relativistic 
Heavy Ion Collider (RHIC) at Brookhaven National Lab- 
oratory is to determine the gluon spin contribution to a 
longitudinally polarized proton (AG), taking advantage 
of the strongly interacting probes available in proton- 
proton collisions [l| . Previous measurements have estab- 
lished the validity of the perturbative Quantum Chromo- 
dynamics (pQCD) description for inclusive mid-rapidity 
7T° [31 and forward 7r° production [3| , and for mid-rapidity 
jet Q and direct photon production [f| , at y/s = 200 
GeV. The double helicity asymmetries for the produc- 
tion of these particles involve gluons in the hard scat- 
tering processes in this pQCD description, and the first 
measurements for n° [1,0] and for jets [H have begun to 
probe AG. 

The RHIC beam polarization and luminosity have sig- 
nificantly improved Q. The statistical uncertainty for a 
double helicity asymmetry measurement is proportional 
to the inverse of P 2 x s/£ for beam polarizations P and 
integrated luminosity C, and decreased by a factor of 5 
from the previously published data from PHENIX [|| 0] • 

In this paper, we first present the cross section for mid- 
rapidity 7T° production for unpolarized proton-proton col- 
lisions at y/s — 200 GeV. These results extend to lower 
and higher px than in previous publications, and we dis- 
cuss an apparent transition region between soft and hard 
scattering; the inclusive cross section is dominated by 
hard scattering, described by pQCD, for px > 2 GeV/c. 
We then present the double helicity asymmetry, All, for 
mid-rapidity 7r° production. We also include measure- 
ments of All at lowpr, below the hard scattering region. 
Finally, our results for px > 2 GeV/c are compared to 
a pQCD calculation that incorporates a model of gluon 
polarization. We present the range that we probe in the 
gluon momentum fraction (x g ) and discuss the constraint 
from these data on AG. 

The PHENIX experiment at RHIC measured 7r°'s via 
71-0 — * 77 decays using a highly segmented (A77 x A<f> ~ 
0.01 x 0.01) electromagnetic calorimeter (EMCal) cov- 
ering a pseudorapidity range of \rj\ < 0.35 and azimuthal 
angle range of Ac/> = ir. The ir° data in this analy- 
sis were collected using two different trigger conditions. 
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A minimum bias (MB) trigger was defined by the co- 
incidence of signals in two beam-beam counters (BBC) 
with full azimuthal coverage located at pseudorapidities 
±(3.0 — 3.9) [lj|. The cross section for events selected by 
the MB trigger was 23.0 mb (about half of c™ ei ) with a 
systematic uncertainty of ±9.7% , derived from vernier 
scan results Q and the variation of MB trigger efficiency 
for subsequent years. Higher px data were collected using 
the coincidence of the MB trig ger and an EMCal-based 
high px photon trigger 0, [Tlf . with efficiency ~ 5% at 
p T (7T°) - 1 GeV/c and ~ 90% for p T (n°) > 3.5 GeV/c. 
The collision vertex was required to be within \z\ < 30 
cm along the beam axis, based on the time difference 
between the two BBC detectors. The ir° acceptance is 
uniform over this interval. The analyzed data sample of 
the 2005 run corresponds to an integrated luminosity of 
2.5 pb- 1 . 

Details of the unpolarized cross section analysis tech- 
nique are described in 0, [Hj]. The background contri- 
bution under the 7r° peak in the two-photon invariant 
mass distribution varied from 80% in the lowest 0.5-0.75 
GeV/c px bin to less than 8% for p T > 4 GeV/c. The 
7T° spectrum was corrected for overlapping decay photon 
showers in the EMCal, based on Monte Carlo simulations 
confirmed with test beam data [l2j]. Below &px(ir°) of 12 
GeV/c the correction is less than 4%, and for px(n°) = 20 
GeV/c the correction is ~25% and ^70%, for two differ- 
ent EMCal subsystems @ . The systematic uncertainty of 
the measurement (excluding the 9.7% uncertainty from 
the MB trigger cross section) varied from ~ 7% at px ~ 1 
GeV/c to ~ 16% for the highest px bin. 

Figure [1] presents the cross section results for mid- 
rapidity 7T° production at y/s = 200 GeV, versus px, 
from p T =0.5 GeV/c to p T =20 GeV [ij]. Points are 
plotted at the average px for each bin. The pQCD 
prediction, at next-to- leading order, is shown for the- 
ory scales [i = px/2, Pt and 2px, where /j, represents 
equal factorization, renormalization, and fragmentation 
scales [ll GJ| . The CTEQ6M parton distribution func- 
tions [l6| and KKP set of fragmentation functions [lTj 
are used. These data extend the published cross section 
data at both low and high px, and are consistent with 
previously published results 0, [ll|. From px—2 GeV/c 
to 20 GeV/c, the NLO pQCD calculation describes the 
data over a change in cross section of seven orders of 
magnitude. 
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The inset to Fig. [T] shows the lower px region in more 
detail including high precision data for the charged pion 
cross section from [18[ ■ The data show a transition in the 
px dependence of the cross section, from exponential to 
a power law dependence, in the region px rsI— 2 GeV/c. 
In order to estimate possible contamination from non- 
perturbative physics in the higher pt data, an exponen- 
tial function (~ e ~ apT ) representing a non-perturbative 
component, is fit to the charged pion spectrum in the 
region px=0. 3 to 0.8 GeV/c (only the lowest pt 7r° data 
point is in this range) and extrapolated to the higher 
Pt region. The exponential fit for the low px region 
gives a = 5.56±0.02 (GeV/c)- 1 , with x 2 /NDF = 6.2/3. 
Only statistical uncertainties for the charged pion data 
were used in the fit. The dominant systematic uncer- 
tainty for the points in the fitted px range is a ~ 12% 
normalization uncertainty (excluding the normalization 
uncertainty from the MB trigger cross section). Beyond 
about pt—^ GeV/c, the data lie above this single expo- 
nential. The fraction of the exponential contribution to 
the data for the 2-2.5 GeV/c px bin is found to be less 
than 10%, with a negligible contribution for higher px- 
This is the basis for applying the pQCD formalism to the 
double helicity asymmetry data with px > 2 GeV/c. 

For the 2005 run, each collider ring of RHIC was filled 
with up to 111 bunches in a 120 bunch pattern, spaced 
106 ns apart, with predetermined patterns of polariza- 
tion signs for the bunches. Spin rotators, sets of four 
helical dipole magnets on each side of PHENIX, rotate 
the polarization orientation from vertical, the stable spin 
direction in the RHIC arcs, to longitudinal at the in- 
teraction point p^j . Beam helicity asymmetries are ob- 
tained by tagging the polarization signs of the bunches 
for each event. The bunches for one beam alternate in 
polarization sign, and pairs of bunches alternate in sign 
for the other beam. In this way data for all combina- 
tions of beam helicity arc collected at the same time, 
and the possibility of false asymmetries due to changing 
detector response versus spin state are greatly reduced. 
Each RHIC fill, typically lasting 8 hours, used one of four 
bunch spin patterns. 

The beam polarizations for 2005 were measured using 
fast carbon target polarimeters (20| , normalized by abso- 
lute polarization measurements made during 2005 by a 
separate polarized atomic hydrogen jet polarimeter [2lj . 
The beam polarizations, from luminosity- weighted 
averages over 104 RHIC fills used in the analysis, were 
(P B )=0.50±0.002(stat)±0.025(systB)±0.015(systG) 
and (P y )=0.49±0.002(stat)±0.025(systY)±0.015(systG), 
for "Blue" and "Yellow" RHIC beams, respectively, for 
the bunches colliding at PHENIX. The systematic 
uncertainties have been separated into uncorrelated 
uncertainties for each beam, "systB" and "systY" , and a 
global systematic uncertainty "systG" , which is common 
for both beams and comes from systematic uncertainty 
in jet polarimeter measurements 22j. For comparison, 
the polarizations in the 2004 run were 0.44±0.08(syst). 

Local polarimeters based on very forward neutron pro- 
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FIG. 1: The neutral pion production cross section at \fs = 
200 GeV as a function of px (squares) and the results of NLO 
pQCD calculations for theory scales fi = pr/2 (dotted line), 
Pt (solid line) and 2pT (dashed line), see text for details; note 
that the error bars are smaller than the points. The inset 
shows, in addition to n°, data for (tv + + tt~)/2 (solid circles), 
and a fit of charged pion data to an exponential function for 
Pt < 0.8 GeV/c (dashed line). The bottom panel shows the 
relative difference between the data and theory for the three 
theory scales. Error bars are quadratic sums of experimental 
statistical and systematic uncertainties (the 9.7% normaliza- 
tion uncertainty is not included). 



duction (production angle 0.3-2.5 mrad) [6|,[23| were used 
to set up and monitor the beam polarization orientation 
at PHENIX. The polarimeters monitor the transverse po- 
larization of each beam at PHENIX, which can be com- 
pared to the beam polarization measured by the RHIC 
polarimeters where the polarization direction is vertical. 
The local polarimeters were calibrated by turning off the 
spin rotators around PHENIX, and measuring the re- 
sponse of the local polarimeters with the beams verti- 
cally polarized. For the longitudinal polarization data, 
the beams showed a measurable transverse polarization, 
with (P T /P) B =0.10±0.02 and (P T /P) Y = 0.14 ±0.02, 
with Px/P referring to the fraction of transverse polar- 
ization of each beam. The polarization directions, as 
determined by the spin rotator settings and as measured 
by the local polarimeters, remained constant over the 
run. The product of the beam polarizations P B ■ P Y 
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is required for the double helicity asymmetry measure- 
ment. The average transverse component of the prod- 
uct was (Eg ■ P£)/(P B ■ P Y ) < (Pt/P) B ■ {Pt/P) Y = 
0.014 ± 0.003; the average of the polarization product 
over the run was (P B ■ P Y ) = 0.24 with a systematic 
uncertainty of ±9.4%. 

The double helicity asymmetry All is the difference 
of cross sections for the same versus opposite beam helic- 
ities, divided by the sum. Experimentally, for inclusive 
7T° production, it can be determined as: 

<=, ' Yl - N++ - R - N+ -; R=*2±, (1) 



5 stat ALL. Therefore, a limit was determined for the Att 
contribution to A LL of 0.07 • 5 stat Aj 



ill- 



where N is the number of 7r°'s measured in PHENIX 
from the colliding bunches with the same (++) and op- 
posite (H — ) helicities, and R is the relative luminosity 
between bunches with the same and opposite helicities. 
Here we neglect the parity-violating difference in cross 

section between (++) < -> f ) and (-1 — ) <-> ( — h) beam 

helicity configurations [24J. All was calculated for each 
fill in order to reduce systematics from variation in beam 
polarizations and in R for different fills. Even and odd 
crossings were handled by separate high px photon trig- 
ger electronics chains. To avoid possible detector bias, 
All was also determined separately for the even and odd 
crossings. Final asymmetries were averaged, and cor- 
rected for the asymmetry of the background under the 
7T° peak in the two-photon mass distribution, as in 

The relative luminosity ratio R is obtained from the 
minimum bias triggers (MB) discussed above. Scalers 
keep track of the number of live triggers for each bunch 
crossing. Single beam background was < 0.05%, as mea- 
sured from non-colliding bunches, and contributes neg- 
ligible systematic uncertainty to the measured R. We 
also measured the double helicity asymmetry of the rel- 
ative luminosity scaler counts, by normalizing using zero 
degree neutral particle production as measured by zero 
degree calorimeters (ZDC) [25| . No asymmetry was ob- 
served. This gave a limit on an asymmetry bias in the 
measurement of (L4J L |bi as < 2 x 10~ 4 , and a limit on 
the systematic uncertainty for the measurement of rela- 
tive luminosity giving SAl L \ R < 2 x 10~ 4 . These limits 
also include the effects from the pileup of two or more 
collisions in a crossing, calculated at < 4% of the cross- 
ings. The BBC and ZDC monitors observe the pileup 
at significantly different rates, and therefore the limits 
above, from comparing BBC and ZDC counts, include 
these uncertainties. 

A transverse double spin asymmetry Att, the trans- 
verse equivalent to Eq. ([1]), can contribute to All 
through the 1.4% transverse component of the product of 
the beam polarizations discussed above. Although Att 
has been postulated to be extremely small, ~ 10~ 4 f26| . 
it has not been previously measured. We measured Att 
in a short run with transverse polarization. Att(pt) was 
consistent with zero within statistical errors [HI; the er- 
rors were 5 times larger than the uncertainties for All, 
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FIG. 2: The double helicity asymmetry for neutral pion pro- 
duction at yfs = 200 GeV as a function of pr (GeV/c). Error 
bars are statistical uncertainties, with the 9.4% scale uncer- 
tainty not shown; other experimental systematic uncertainties 
are negligible. Four GRSV theoretical calculations based on 
NLO pQCD are also shown for comparison with the data (see 
text for details.) 



Figure [2] presents the measured double helicity asym- 
metry in 7T° production 13]. A scale uncertainty of 9.4% 
in A\ h due to the uncertainty in beam polarization is not 
shown. The other systematic uncertainties are negligible, 
as discussed above, and checked using a bunch polariza- 
tion sign randomization technique, and by varying the 7T° 
identification criteria @. Data for pt > 1 GeV/c were 
obtained from the high pt photon triggered sample. For 
Pt below 1 GeV/c, due to low efficiency for the high pt 
photon trigger, we used the MB data sample. In the low 
Pt region, where the cross section shows an exponential 
behavior, the helicity asymmetry is A£ L =0.002±0.002, 
for the data in the range pt = 0.5 — 1 GeV/c. For the 
higher pt region, the four curves in Fig. [5] show calcula- 
tions of A\ L , using NLO pQCD with fi = pt{tt°), that 
reflect the range of gluon polarizations allowed by inclu- 
sive deep inelastic scattering (DIS) data. The calcula- 
tions are based on the GRSV model, where "std" was 
the best fit to inclusive DIS data [27j | . For momentum 
fraction x, AG(x) — G + (x) — G~(x) refers to the gluon 
helicity distribution, and G + (x) and G~(x) refer to the 
gluon densities for + and — helicities in a + helicity pro- 
ton. The first moment of the gluon helicity distribution, 
l 

J AG(x)dx, for the "std" parameterization is AG=0.4, 
o 

at the scale Q 2 =l GeV 2 . The other three curves are cal- 



fi 



culations based on this best fit, but use at Q 2 = 0.4 GeV 2 
the function AG(x) — G(x),0, —G(x), where G(x) is the 
unpolarized gluon distribution. The gluon distribution 
at the input scale is evolved to the scale Q 2 = p^-ir ). 

In order to explore the impact of the new data on the 
sensitivity to the polarized gluon distribution, we have 
compared the data with a set of All{pt) curves corre- 
sponding to different AG(x) between AG(x) — —G{x) 
and AG(x) = G(x) at Q 2 = 0.4 GeV 2 . We used the 
data for px > 2 GeV/c, which appear to have little con- 
tamination from soft physics as discussed earlier. The 
most likely Xn for PHENIX tt° data in each px point is 
~ xt/0-8 |3(j, where xt — pt/(\/s/2). For the mea- 
sured pr range 2-9 GeV/c, the range of x g in each bin is 
broad, and spans the range x g = 0.02 — 0.3, as calculated 
by NLO pQCD [p]. 



Figure [3] shows the 
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X versus 
an integral 

of AG over the probed x g range. Only experimental 
statistical uncertainties are used to calculate % 2 , and no 
theoretical uncertainties are included. It is important 
to note, that although the range of the first moment 
explored represents ~ 60% of the full integral, this 
reflects using a specific model for the gluon polarization. 
For example, a gluon polarization model with a crossover 
from positive to negative gluon polarization within our 
x g range would yield a small average asymmetry for each 
point. Also, other models can generate larger or smaller 
contributions from the gluon spin in the unmeasured 
region of x g . 

These data are sensitive to the first moment of the 
polarized gluon distribution. Using the GRSV model, 
we find that the gluon polarization contribution to the 
proton spin (1/2) in the probed x g range is constrained 
between —0.9 and +0.5, for x 2 — Xmin=9> representing 
a "3er" limit (a "ltr" limit would give a constraint be- 
tween 0.07 and 0.3). The extremes of gluon polarization 
are ruled out, modulo the above remarks, with the con- 
fidence level for "AG = ±G" of less than 10 -6 . Large 
positive gluon polarization [28| was proposed shortly af- 
ter the discovery that the quark contribution to the pro- 
ton spin was small [2^ |. with the suggestion that such a 
large gluon polarization would mask a "bare" quark po- 
larization. For "std" and "AG = 0" , the confidence lev- 
els are 20-21% and 12-13%, respectively, for the range 
of ±9.4% scale uncertainty of the measurement. Semi- 
inclusive DIS measurements (32j have also presented 
data on AG in a limited x g range and its comparison 
with various AG scenarios. 

The two minima in Fig. [3] reflect the quadratic contri- 
bution of the gluon polarization to All, from the gluon- 
gluon scattering subprocess for ir° production. The sym- 
metry between the two minima is broken by the quark- 
gluon scattering subprocess, where the gluon polarization 
contributes linearly to All- The quark-gluon subprocess 
is emphasized at higher p?, which will become accessible 
with additional running at high polarization and lumi- 
nosity. 
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FIG. 3: The distribution of the measured data plotted 
versus the value of the first moment of the polarized gluon 
distribution (solid line) in the x g range from 0.02 to 0.3 cor- 
responding to our 7T° data in pr bins from 2 to 9 GeV/c. 
Dashed and dotted lines correspond to —9.4% and +9.4% 
variation in All normalization related to the beam polar- 
ization uncertainty, the dominant systematical uncertainty of 
our data. Only statistical uncertainties were used for each 
curve. Arrows indicate AG corresponding to the different 
polarized gluon distributions discussed in the text. 



To summarize, we have presented the unpolarized 
cross section and double hclicity asymmetries for ir° pro- 
duction at mid-rapidity, for proton-proton collisions at 
■s/s = 200 GeV. We observe an apparent transition re- 
gion in the cross section, for pt ~1 to 2 GeV/c, with 
the cross section described by an exponential in pt be- 
low pr ~ 1 GeV/c, and with the cross section described 
by the pQCD prediction for px — 2 to 20 GeV/c, over 
seven orders of magnitude in cross section. The results 
for All in the pQCD region, which we take as pr > 2 
GeV/c, constrain the gluon polarization in the proton 
significantly. The range probed is x g = 0.02 to 0.3, for 
the gluon momentum fraction. Using one representative 
model for the gluon polarization, GRSV |27j, which as- 
sumes no crossover in gluon polarization versus x g , we 
present a map of x 2 versus the first moment of the polar- 
ized gluon distribution in the measured region. From this 
study, the present data rule out extreme values of gluon 
polarization suggested after the surprise of the EMC re- 
sult that the quarks (and anti-quarks) contribute little 
to the spin of the proton (29|, but allow significant con- 
tribution from the gluon spin to the proton spin. 
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